Optimal generation of indistinguishable photons 
from non-identical artificial molecules 
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We show theoretically that nearly indistinguishable photons can be generated with non-identical 
semiconductor-based sources. The use of virtual Raman transitions and the optimization of the 
external driving fields increases the tolerance to spectral inhomogeneity to the meV energy range. 
A trade-off emerges between photon indistinguishability and efficiency in the photon-generation 
process. Linear (quadratic) dependence of the coincidence probability within the Hong-Ou-Mandel 
setup is found with respect to the dephasing (relaxation) rate in the semiconductor sources. 
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Single-photon sources (SPSs) are fundamental devices 
in quantum communication pQ and linear-optics quantum 
computation [2]. Essentially, a SPS consists of an atomic- 
like system that can be deterministically excited and thus 
triggered to emit single-photon wavepackets into a prefer- 
ential mode. In addition, all photons have to be emitted 
in the same quantum state, in order to maximize the 
visibility of two-photon interference. 

So far, single-photon generation has been demon- 
strated both with atomic HH03IB], molecular 0151 US, 
and solid-state systems [TOl HU EEU HI HI EES] • Amongst 
the latter, self-assembled semiconductor quantum dots 
(SAQDs) |16j are particularly promising, for they can be 
controllably coupled to optical microcavities [T7], so as 
to greatly enhance both the photon emission rate and 




FIG. 1: (Color online) Emission scheme of indistinguishable 
photons from non-identical SPSs. Each source (A/B) consists 
of an AM (gray), doped with an excess electron (red/blue), 
and coupled to an optical MC (black). The effect on the 
photon wavepackets (green) of the differences between A and 
B, in terms of energies and oscillator strengths of the optical 
transitions, are compensated by properly tailoring the excit- 
ing laser pulses (straight arrows). 



the collection efficiency. As a major downside, SAQDs 
come with a finite dispersion in terms of size, shape, 
and composition. This typically gives rise to spectral 
inhomogeneities in the meV range, i.e., orders of magnti- 
tude larger than the homogeneous linewidths (few /ieVs) 
of the lowest exciton transitions at cryogenic tempera- 
tures. Photons spontaneously emitted by two distinct 
SAQDs therefore tend to be completely distinguishable: 
this could impede scalability and ultimately limit the po- 
tential of semiconductor-based SPSs. 

In spite of its relevance to any solid-state approach, 
the problem of generating indistinguishable photons with 
non-identical emitters has still received limited attention. 
In this Letter, we assess the possibility of compensating 
the effects of such differences between two SPSs by suit- 
ably tuning the exciting laser pulses. To this aim, we 
combine a density-matrix approach |18j - to simulate the 
system dynamics and compute the coherence functions of 
the emitted radiation - with a genetic algorithm |25j - to 
optimize the external driving fields. As a crucial point, 
the photon generation results from a (virtual) Raman 
transition. Raman transitions have already been pro- 
posed as means to avoid the classical uncertainty on the 
initial time of the emission process (the so-called time- 
jitter) [T5J [20] and to tune both the temporal profile and 
the central frequency of the emitted wavepacket [5] . Here 
we show that such flexibility can be exploited to gener- 
ate two nearly indistinguishable photons, in spite of the 
spectral differences between their respective semiconduc- 
tor sources. 

We specifically consider the case where each source is 
represented by two coherently coupled SAQDs (often re- 
ferred to as artificial molecule, AM) , embedded in an op- 
tical microcavity (MC) (Fig. [TJ. The AM is doped with 
an excess electron, whose levels 1 and 2 [Fig. |2ja)] define 
a pseudo-spin, corresponding to the hybridized - bond- 
ing and antibonding - states of the two dots [211 f22] : 
together with the lowest charged exciton states (3 and 
4), these form a double A— scheme [53]. There, a Raman 
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FIG. 2: (Color online) (a) Relevant level scheme for each 
negatively-charged AM. The levels 1 and 2 correspond to the 
bonding and antibonding states of the excess electron; these 
are optically coupled to the two lowest charged exciton states 
3 and 4. The off-resonant laser pulse (<jJl) induces a virtual 
Raman transition from 1 to 2; this process results in the cre- 
ation and emission of a photon from the fundamental mode of 
the MC (iv c ). (b) Schematics of the Hong-Ou-Mandel setup. 
The photons are emitted by the two sources A and B in the 
input modes of a balanced beam-splitter (BS), whose output 
modes are coupled to the photodetectors C and D. 

transition can be induced by an off-resonant laser pulse 
(u>l), that drives the AM from 1 to 2, while generating a 
cavity photon (w c ). 

We simulate the dynamics of the two sources (A and B) 
within a density-matrix approach [18]. The state of each 
AM-MC system is denoted by \j,n), where j = 1, ... ,4 
specifies the AM eigenstate, whereas n represents the oc- 
cupation number of the cavity mode. The time evolution 
of the density matrix within such Hilbert space is given 
by p — i[p,Hmc + Hi] + Cp (rotating- wave approxima- 
tion, H = 1). The coupling of the AM with the MC is 
accounted by 

H M c= ^2 E 9jk(<7jk^ + acr] fc ), (1) 
3=1,2 k=3A 

where a is the annihilation operator of the cavity photon 
and Ojk = \j){k\ the ladder operator acting on the AM 
state. The coupling of the AM with the exciting laser is 
given by 

1 Nl 

51 = 2^1 E WCe-"'*'^ + H.c), (2) 

p=l j=l,2 fc=3,4 

where ^ k (t)/gjk = ex P{(* — ^o) 2 /2 (T p} gives the time 
envelope of the p— th laser pulse, 5^ k = e& — 6j — u L , and 



£j the AM levels. The interaction of the AM-MC system 
with the environment is given by the superoperator 

C = Y, L (V^ P j)+T, L (^ik^k)+HV^a) 1 (3) 

3 j,k 

where L(A)p = (2ApA^ - A^Ap - pA*A)/2 and Pj = 
\j){j\. The three terms in Eq. [3] account for pure de- 
phasing (with rate 7), radiative (jk — 13, 14,23,24) and 
non-radiative (jk = 12, 34) relaxation of the AM, and 
cavity-photon emission, respectively. In the following, 
we assume for simplicity that Tjk = T r for all the radia- 
tive relaxation processes of the AM, and Tjk = T nr for 
the non-radiative ones. 

The degree of indistinguishability between the photons 
generated by the two sources can be measured within 
the Hong-Ou-Mandel setup [24 [Fig. (2j»]. There, if 
two indistinguishable photons enter the input ports (A 
and B) of a balanced beam splitter, two-photon inter- 
ference results in a vanishing probability of a coincidence 
event (Pcd = 0) in the two photodetectors at the output 
modes (C and D). For distinguishable photons, instead, 
Pcd — 1/2. The coincidence probability can thus be re- 
garded as a measure of the photon indistinguishability. 
Formally, Pc d can be connected to the dynamics of the 
A and B sources through [T^] Pcd = F\ — F%, where 
F p =x,2 — k 2 J dt J dr f p (t,r) and 

/i(i,r) = [n A (t + T)n B (t) + n B (t)n A (t + T)}/4, 
f 2 (t,r) = Re{[G^(t,t+T)rG^(t,t + r)}/2. 

Here, G^\t,t + r) = (a\,(t)a x (t + t)) and n x (t) = 
(oJ.(t)a x (t)) are, respectively, the first-order coherence 
function and cavity-mode occupation corresponding to 
the sources x — A, B. 

In order to maximize the overlap between the 
wavepackets of the photons emitted by A and B, we 
optimize the driving laser pulses and the frequency of 
the cavity mode, i.e., the vector X = (CIq , t$ , a n , ui c ). 
For two given sources [each characterized by the vector 
Y = (£j,gjk,~/jiFjk), with Y,4 ^ Yb], the suitability of 
each set of laser pulses for the generation of two indis- 
tinguishable photons is quantified by a fitness function 
F(X A ,X B \Y A ,Y B ) > 0, defined as 

F=(F 2 /F 1 )g(P e A ,P e B ). (4) 

Here, < F2/F1 < 1 measures the degree of indistin- 
guishability between the two photons, whereas g accounts 
for the statistics of the two SPSs; more specifically, it 
imposes a penalty on the vectors X x corresponding to a 
photon-emission probability P e x = J n x (t)dt far from 1. 
We identify the best solution (X°) with the vectors that 
correspond to the maximum value of the fitness func- 
tion, for the given values of the physical parameters that 
characterize the SPSs (Y x ): T M = ^(X^X^Y^, Y B ). 
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FIG. 3: (Color online) Maximized fitness function Tm as 
a function of the differences between the A and B sources 
in terms of frequencies (Aab = 1^32 — ^32) an( A oscilla- 
tor strengths (gAB = gfk/dfk) °f the optical transitions. 

A/ B 3 1 

The remaining parameters are: r nr = T r ' = 10 ps , 
y A/B = 10 -a VB -i t k a/b = io-ipB-i, gf k = .120meV. 
Having found no evidence of an increased indistinguishabil- 
ity arising from non-gaussian laser-pulse profiles, we set here 
N L = 1. 

Hereafter, we take: g = 8(Pf - P l e ) ■ 0(P e B - P*), with 
8 the Heaviside function. Therefore, provided that both 
P e A and P e B are larger than the threshold P* , 1 — T co- 
incides with the coincidence probability Pcd, normal- 
ized to the joint emission probability of the two sources 
Pab — P£ • Pe — Pi- To this aim, we combine the 
density matrix approach with a genetic algorithm |25j : 
this allows to efficiently explore, within a large parame- 
ter space, the particularly complex landscape induced by 
the penalty function g. 

The distinguishability between the photons emitted 
by the two sources mainly arises from the differences 
between their transition energies (^>fk ^ uif k , being 
ujjk = Cj — €k) and oscillator strengths (gf k ^ gfk)- ^ n 
Fig. [3] we show the dependence of the maximum fitness 
for optimized laser pulses, Tm, on 9ab = gfk/dfk an< ^ 
Aab — ^32 — ^32 [26 , for a threshold emission prob- 
ability P l e = 0.9. The fact that T M < 1 even for the 
case of identical sources (gAB = 1 an d Aab = 0) is 
due to the presence of dephasing and non-radiative re- 
laxation (see below). By adjusting the laser pulses, mod- 
erate differences between the oscillator strengths of the 
two sources (gAB ^ 0.9) can be efficiently compensated, 
whereas decreases of Tm arising from mismatches be- 
tween the optical-transition energies of A and B in the 
meV range can be limited to a few percent. 

A couple of general comments are in order. On the one 
hand, stimulated Raman processes do increase the toler- 
ance to inhomogeneities between the SPSs from the /xeV 
- as is the case with spontaneous emission - to the meV 
range. On the other hand, such inhomogeneities cannot 
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FIG. 4: (Color online) (a) Dependence of the maximum fit- 
ness function Tm on the dephasing rate 7 (upper axis, green 
squares) and on the non-radiative relaxation rate T2i (lower 
axis, red triangles), for g AB = 1 and A^a = 1-25 meV. The 
green (red) dotted line represents the best fit of Tm linear 
(quadratic) in 7 (r2i). (b) Tm as a function of the thresh- 
old P* in the photon-emission probability for gAB = 1, with 
Aab ~ and 1.5 meV. 



be completely compensated by properly tuning the fre- 
quencies of the driving laser pulses, as could be naively 
expected on the basis of a simple energy-conservation re- 
lation. These limitations are partially due to the effects 
of decoherence; besides, the requirement that the two 
sources emit photons with a high probability seems to 
conflict with that of maximizing their indistinguishabil- 
ity In the following we investigate separately these two 
aspects. 

The coupling between the carriers confined in each AM 
and the phonons gives rise to energy relaxation and de- 
phasing. These two contributions are considered sepa- 
rately in Fig. Qa) , where we report Tm &s a function of 
7 (with = 0, squares), and of T21 (with 7 = 0, trian- 
gles), within realistic ranges of values for these rates. As 
in the case of two photons sequentially emitted by the 
same source |19j . the coincidence probability increases 
linearly with 7 (dotted green curve). A stronger depen- 
dence is found with respect to the non-radiative relax- 
ation between the two lowest states in the A— scheme (a 
fit by a quadratic function of T2i is given by the dotted 
red line). We note, however, that the value of T21 can be 
strongly reduced by suitably engineering the AM geom- 
etry and by the application of external fields [S7J HE] • 

In order to investigate the relation between the emis- 
sion probability of the A and B photons and their in- 
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FIG. 5: (Color online) (a) Dependence of F 2 /Fi = 1 - P CD 
on the shift of the cavity and laser frequencies from their op- 
timized values, Awf fl = Ua^ B — u>a^ B '° , for Aab = (red) 
and Aab = 1.5 me V (gray), with qab = 0.88. Solid, dashed 
and dotted lines correspond respectively to: a = c; a = L, 
[3 = A; a — L, P — B. (b) Photon-emission probabilities, 
P^ (solid lines) and Pf (dotted), as a function of Auj c , for 
Aab = (red) and Aab = 1.5 meV (gray), with gAB = 0.88. 



distinguishability, we have computed Tm as a function 
of P*. In Fig. |4^b) we report two representative cases: 
Aab = (orange squares) and Aab = 1.5 meV (gray 
triangles), with gAB = 1- The photon emission proba- 
bility and the degree of indistinguishability are clearly 
anti-correlated, thus demonstrating the existence of a 
trade-off between the two requirements. The dependence 
of Tm ° n Pe> though evident even for identical sources 
(Aab = 0), is more pronounced in the presence of a spec- 
tral mismatch (Aab — 1.5). The above curves are rela- 
tively insensitive to differences in the oscillator strengths 
(gAB ^ 0.9, not shown here). 

We finally comment on the robustness of the solution 
with respect to small departures of the control parame- 
ters from their optimal values. As a representative ex- 
ample, we report the dependence of F%/Fi — 1 — Pcd 
on the cavity frequency uo c = uj^ — ojf , referred to its 
optimized value oj° [Fig. [5] (a)]. The robustness of the 
photon indistinguishability with respect to non-optimal 
parameters decreases for incresing spectral mismatches 
between the two sources. Besides, our simulations show 
a stronger dependence of Pcd on the laser (solid and dot- 
ted lines) than on the cavity frequencies (dashed). This 
feature can be traced to the fact that the uncertainty on 
the cavity frequency (~ 1/k = 10 ps) is larger than that 
on the laser frequency (typical duration of the optimized 
laser pulses are few tens of ps). As to P^ B [Fig. [5] (b)] , 
we note that the value corresponding to the optimized 
parameter (Au> c = 0) is close to the mimimum allowed 
value (Pi): this provides a further indication on the ex- 
istence of a trade-off between photon indistinguishability 
and emission probability. Besides, the existence of a spec- 
tral mismatch (gray lines) makes also the fulfilment of the 
requirement P^,Pf > P* more sensitive to the tuning 
of the physical parameters. The required precision is of 
the order of a few /ieV for Aab = 1.5 meV. 

In conclusion, we have investigated the degree of indis- 



tinguishability between single photons emitted by non- 
identical artificial molecules. We found that the use of 
virtual Raman transitions, combined with the optimal 
design of the driving laser pulses, increases the tolerance 
with respect to the spectral inhomogeneities between the 
sources up to the meV energy range. Defined power-law 
dependences of the coincidence probability on the de- 
phasing and (non-radiative) relaxation rates are identi- 
fied. Besides, unlike the case of identical sources, a trade- 
off emerges between the requirements of maximizing the 
emission efficiency and the photon indistinguishability. It 
is finally worth mentioning that the use of a pseudo-spin 
within the A-level scheme offers promising possibilities 
|31j . These include the use of quantum-confined Stark 
effect [2pJ for triggering the photon emission [3D], and 
the entangling of carrier spins localized in remote (and 
thus non-identical) semiconductor systems [32] . The lat- 
ter goal requires the combination of a spin-photon entan- 
gling [33] process with two-photon interference within the 
Hong-Ou-Mandel setup. 
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